Low energy valence photoemission in the Anderson impurity model for Ce compounds 



0\ 



3 



o 

CZ3 



J.D. Lee 

Max-Planck Institut fur Festkorperforschung, 
Heisenbergstrasse 1, D-70569 Stuttgart, Germany 
(February 7, 2008) 

The valence level photoemission spectra in the Anderson impurity model for Ce compounds at 
zero temperature are studied as a function of the photon energy u). Most of former studies on Ce 
compounds are based on the sudden approximation, which is valid in high energy region. For the 
photoemission in the adiabatic limit of low energy region, one should consider the dipole matrix 
elements and the hole-induced photoelectron scattering potential. We can manage it by combining 
the time-evolution formalism and the 1 /Nf scheme in a large /-level degeneracy Nf . This gives the 
exact results as Nf — » oo. In view of experiments on the valence photoemission, two contributions 
of 4/- and band emissions are mixed. We study the separate 4/ and band contributions (from 
Ce 5d) and total emission including the interference between two on an equal footing with varying 
the photon energy. In the 4/-emission case, we also explore the effects of hole-induced scattering 
potential of the photoelectron with respect to u). Its effects are found very similar to the core 
level photoemission in shake down case with a localized charge transfer excitation. Additionally, 
we examine the adiabatic-sudden transition in valence level photoemission for the present localized 
system through the simplified two-level model. 
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I. INTRODUCTION 



The Anderson impurity model (AIM) was originally 
proposed to discuss the property of magnetic impuri- 
ties in nonmagnetic metalsa. After that, AIM has been 
widely applied to the analysis of spectroscopic data-,for 
/ and d electron systems, i.e. rare earth compoundsa or 
transition metal compounds^, where electron states are 
treated to be an impurity and they are hybridized with 
the valence or conduction electron states. Also, AIM has 
been often used to describe the Ce mixed-valence com- 
pounds, where one considers the / level on one atom 
and its interaction with the conduction bands. In inves- 
tigations of Ce and its compounds, the basic question 
concerns the nature of a 4/ electron and other electronic 
states and how they mix with the 4/ state. Much of the 
interests are therefore imposed on the properties of the 
4/ states, i.e. the occupation, position, width, coupling 
to the metallic band, intra-atomic Coulomb interaction 
and so on. There were numerous studies of thermody- 
namic and transport properties for them, which has been 
followed by the electron-spectroscopy studies!! 

The photoemission spectroscopy (PES) is a very useful 
tool for studying the electronic structure of matters and 
could have provided a lot of insights also for Ce stud- 
ies. But it's worth noting PES cannot always give a 
simple answer about the underlying electronic structure 
because the photoelectron may perturb the system left 
behind. Actual description of theoretical PES is quite 
complicated and therefore the sudden approximation is 
frequently used, where the photoelectron is assumed de- 
coupled from the remaining solido. The sudden approx- 
imation becomes exact when the kinetic energy of the 
emitted electron gets large infinitely. 

Gunnarsson and Schonhammem have extensively stud- 
ied the electron spectroscopies for Ce compounds, 



i.e. core level photoemission, x-ray absorption, and 
bremsstrahlung isochromat spectroscopy as well as the 
valence photoemission. For the information of the po- 
sition and width of /-level in the compounds, the va- 
lence photoemission has often been used. Through their 
studies of valence photoemission in /-emission channel, 
they reproduced two-peak structure in Ce compound us- 
ing AIM consistent with the experiments. In the earlier 
evolution stage for Ce materials, it was found that the 
valence PES shows just a gfflgle /-related structure 2- 
3 eV below the Fermi leveDtl. Later PES experiments 
have demonstrated the 4/ spectrum has the -additional 
structure interestingly near the Fermi leveEilj. Subse- 
quently, it was shown that the particular structure is due 
to the Kondo resonance singlet characterized by the small 
energy Tr- . In the actual experiments on valence photoe- 
mission, two contributions of 4/-emission and band emis- 
sion (from Ce 5d or other bands) are mixed. The identi- 
fication of 4/-emission from the experiments is a highly 
nontrivial work. Wieliczka et alJij have done the com- 
parison of spectra taken at the different photon energies 
showing the resonance of 4/-emissionL l 3 and reported the 
additional peak near the Fermi level. Another possibil- 
ity is to assume the behaviors of 4/- and band emissions 
with, respect to the photon energies, especially in 20 — 80 
eVt-X Nevertheless, they could have said nothing about 
the interference between two. These works can motivate 
to explore a more explicit analysis for the interference 
effects of two emission channels with the photon energy 
varied. Gunnarsson and SchonhammeitJ have also stud- 
ied the band emission contribution as well as 4/-emission 
and discussed the interference effects of two. However, all 
their works were within the sudden approximation. 

We consider both the contributions of 4/- and 5d- 
emission on an equal footing by introducing two dipole 
matrix elements, Af(E) and Ad(e,E), where E is the 
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kinetic energy of photoelectrons. In the sudden approxi- 
mation, the dipole matrix elements are normally treated 
constant with E. But in the low energy PES, it can 
be crucial. Combining the time-dependent formalism 
and 1/Nf idea, we can calculate the PES exactly up to 
O(jjj) as the photon energy varies. Nf = oo can be 
a good approximation to Nf = 14 in Ce compounds. 
Then we study the separate contributions of 4/- and 5c?- 
emission and more interestingly the nontrivial interfer- 
ence effects between two. The relative sign or strength 
of Af(E) and A d (e,E) are important. Difference in the 
energy scale of Af(E) and A d (e, E) makes the spectra 
from each channel separate with respect to u> in the low 
energy PES. It is also found that, because the interfer- 
ence contribution has a peak near the Fermi level, the 
4/-derived peak near the Fermi level may be enhanced 
or suppressed in a total spectra. 

It is recently reported that the adiabatic-sudden transi- 
tion due to the photoelectron scattering potential will be 
governed by the characteristic of relevant excitations to 
which the emitted electron couples in the system. When 
the photoelectron couples to the extended excitation like 
piasmon, the sudden transition occurs in very large ki- 
netic energies (~ keV)t3, while to the localized excita- 
tions, the sudden approximation occurs much quickerEH. 
In our study, the extrinsic effects beyond the sudden ap- 
proximation will be considered only through the hole- 
induced scattering potential. Effects like surface or sev- 
eral damping mechanisms will not be taken into account. 
The system of Ce compounds has also included the local- 
ized excitations represented by /°, f 1 , f 2 created from 
the hole potential. We can consider the scattering po- 
tential in the valence PES due to /-hole in the impu- 
rity model within the formalism. The scattering effects 
will not be important for band emission channel. The 
crossover from adiabatic to sudden limit can be also re- 
examined in this localized system, for which we simplify 
the model to have only two relevant levels. The same 
criterion for the transition is found as in the previous 
workO, the energy scale of E — l/(2R 2 ), where R is a 
scattering potential range. 

We organize the paper as follows. Our model and the 
formalism for calculation are given in Sec. II. The sim- 
ple sudden approximation results for separate and both 
channels are described in Sec. III. In Sec. IV, we present 
the model for necessary dipole matrix elements and cal- 
culate the spectra for the separate contributions for 4/- 
and 5<i-emission and for both with respect to the photon 
energies. We also discuss the interference contribution 
between two. In Sec. V, for 4/-emission, we include the 
photoelectron scattering potential within the same for- 
malism. In Sec. VI, we try to reexamine the adiabatic- 
sudden crossover in the system by way of the simplified 
two-level model. In See. VII, we give the discussion and 
conclusion. 



II. MODEL AND FORMALISM 

As mentioned in the introduction, we consider the AIM 
Hamiltonian in the energy basis used in GunmLcason and 
Schonhammer's discussion for Ce compoundsotJ, 

H = Ho + A, 

Ho = I E ^ E ^EudE + I £ ^l^ de + e f J2 n » 

where i\) Ev {^Ev) is a photoelectron operator, e denotes 
the 5d conduction states, e/ describes the impurity 4/ 
level, and V(e) is a hybridization matrix element between 
the conduction states and localized / level. | V(e) \ 2 can be 
modelized to have a semielliptical form symmetric with 
respect to t F = 0, 7r|V(e)| 2 = 2V 2 (B 2 - e 2 ) 1 / 2 / B 2 , where 
2B is the bandwidth. A in the Hamiltonian is the dipole 
term describing the photon-matter interaction. The ona 
particle basis used in Eq.(|l]) is introduced by assumingEEl 

5>k* m VW = \V(e)\ 2 6 mm ,, 



1>L = V{e) 



k 



and so v is the orbital and spin magnetic quantum num- 
ber and from v = 1 to v = Nf if we assume the magnetic 
degeneracy Nf of /-level. In Ce, Nf is normally taken as 
14. To apply 1 /Nf idea, we need one subsidiary condition 

1 /2 

that Nf F(e) should be independent of Nf. 

In A, we will generally have two interaction terms due 
to 4/ level and 5d conduction bands, so 

A = J2jdE [A f (E)^ v + A){E)^ Ev ] (2) 



— J2 / dedE [A d (e, E)^ Ev ^ ev + A* d (e, E)^ Eu 



N 



By giving the explicit time-dependency /(r) in A and re- 
defining the dipole matrix elements, Af(E) and Ad(e, E), 



A f (E) -> M f Af(E)f(r), 
A d (e,E)^M d A d (e,E)f( T ), 



i), v>o 



(3) 



(4) 



we use a time-dependent formulation and solve the 
Schrodinger equation for the total Hamiltonian TL. 
We first introduce a state |0), 



N f 



i°)=n n ^iivac), 



(5) 



v £<£f 



2 



where all the conduction electron states below Fermi en- 
ergy are occupied and the / level is empty. For the sim- 
plicity, we keep only the lowest order terms of 1 /Nf be- 
fore and after the photocmission, which means the results 
will be exact as Nf — > oo. 



\E,e,e') 



N 



(6) 



1^^ <2( ^_ 1)1/2 E^^^10), (7) 



\E,e) 



N 



1/2 



(8) 



JVj /2 (JV/-l) 1/a ^ 



E ^L^Vw^|0>, (9) 



let the system evolve for a time of the order 1 jr\ . Then 
we can show the sol ut i on identical to the more conven- 
tional photoemissionllZHIa. r\ is a small positive number 
and gives a life-time broadening in the spectra. In the 
actual calculation, r\ is taken as 0.3 eV (0.01 au). The 
photocmission spectra will now be proportional to 

1(E) = J \d(E,e- 7 T)\ 2 de + J \e(E,e,e';T)\ 2 dede' 

+ J \f(E,e,e',e";r)\ 2 dede'e", (14) 

and we see, due to Mf — > and Md — > 0, 

1(E) = a(E)M 2 + f3(E)Ml + j(E)MfM d , (15) 

where a(E), (3(E), and j(E) correspond to 4/-, 5d- 
emission, and interference between those, respectively. 

III. SUDDEN APPROXIMATION 



\E,e,e',e") = 



N) /2 (N f - iy/ 2 (N f - 2)V2 

X J2 4^4"VVV»^eV'^M0>. (10) 

V^l/'^li" 

Within the above basis set, after time r, the wave func- 
tion |\&(t)) of the system is given by 



|¥(r)> = o(r)|0) + / b(e;r)\e)de + J c(e, e'; r)|e, e')dede' 

+ J d(E,e;T)\E,e)dEde 

+ J e(E,e,e';r)\E,e,e')dEdede' 

+ J f(E,e,e',e";T)\E,e,e',e")dEdede'de". (11) 

The coefficients of |^(r)) can be determined by the time 
dependent Schrodinger equation, 



iA|*(r))=W|*(r)), 

OT 



(12) 



where the initial condition of the state should be corre- 
sponding to the ground state before the photoemission, 
|¥(r = Q)) = |*o>, 



|*(0)) = o(0)|0> + / b(e;0)\e)de + J c(e, e'; 0)\e, e')dede' 

(13) 

and the equations for a(0), b(e; 0), and c(e, e'; 0) are found 
in Ref.15. The coefficients Mf and Md represent the ex- 
ternal field strength. In the present formalism, we solve 
the equation in the limit of Mf — > and Md — ► and 



In the sudden approximation, we normally neglect in- 
dependency of the dipole matrix elements, i.e. Af(E) = 
Mf and Ad(e,E) — M d A' d (e), where £ is a kinetic ep. 
ergy of photoelectron. Gunnarsson and SchonhammerEj 
assumed A' d (e) have the same shape as V(e) for the con- 
duction band emission, which we'll simply follow. 

The AIM has often been studied in the limit of U — 00, 
where it becomes so simple as to allow the analytic solu- 
tions. In our formalism, to neglect |e, e') and \E,e,e',e") 
corresponds to the limit. However the assumption U = 
00 is not really justified because U is just about 5^6 
cVcJ and thus f° and f 2 configurations are energeti- 
cally comparable, i.e. e/ is about —2 ~ —3 eV0 and 
2ef + U » 0. In the calculations, we have always taken 
U = 5.0 eV and e f = -2.5 cV to be 2e/ + U = 0. In 
Fig.jl], we give the simple valence PES results for 4/- and 
5d-emission, respectively. 
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FIG. 1. In the upper panel, the /-derived valence PES is 
provided (Md = 0) and in the lower panel, the conduction 
band emission is given (Mf = 0) for U = 5.0 eV, e/ = —2.5 
eV, V = 0.5 eV, and B = 6 eV. Both calculations are based 
on the sudden approximation. Spectral curves are normalized 
to have the same area. 

It is seen in Fig.[l] that we jaicely reproduce the well- 
known sudden 4/-PES rcsultslla having the double-peak 
structure in the upper panel and also get 5d-PES simu- 
lating the broad structureless conduction band. The 4/- 
PES is especially interesting because of its ample physics. 
The peak well below the Fermi level corresponds to 4/ 
ionization peak, 4/ 1 — > 4/°, and the peak near the Fermi 
level (also called Kondo resonance peak) arise from a 4/ 
hole screened by a 4/ electron (making a bd hole near 
the Fermi level), 4/ 1 — > if 1 . 

In the separate calculations of emission channels, the 
absolute values or signs of Mf and M d are surely irrel- 
evant to the results. There can be however some sub- 
tleties when we consider both emission channels. The 
PES curves drastically change with respect to the rela- 
tive sign or relative ratio of Mf and Md- For the rela- 
tive strength, we parametrize the ratio of \MdA'(0)/Mf\ 
(note \MdA'(0)/Mf\ > 1 does not always mean the band 
emission is dominant over the /-emission). The relative 
sign is related to whether the interference will be con- 
structive or destructive. 




ENERGY (eV) 

FIG. 2. The valence PES including both contributions of 
4/- and 5d-emission is shown. The solid line is the total 
spectra corresponding to constructive or destructive interfer- 
ence, where interference effects (labeled "intf) is added or 
subtracted. In (a), we use |M d A'(0)/M/| = 2 and in (b), 
|Af<jA'(0)/M/ 1 = 6. The used parameters are same in Fig.l. 



As shown in Fig.g, the relative ratio and signs of two 
dipole components are crucial in valence PES. Interest- 
ingly, the interference contributions show a peak near the 
Fermi level, which may enhance the Kondo resonance 
peak from 4/-emission in the constructive case or sup- 
press in the destructive case. In FigJ|(b), even if we can 
see a clear ionization peak near e/, we see only the shoul- 
der structure not a peak near the Fermi level due to a 
strong destructive interference. These destructive inter- 
ferences may be one of the reasons that several groups in 
the former stage have failed to see a double-peak struc- 
ture and reported only an ionization peak near e/. 

IV. LOW ENERGY VALENCE PHOTOEMISSION: 
EFFECTS OF DIPOLE MATRIX 

In the last section, we have illustrated the sudden ap- 
proximation results valid in the high energy PES. To see 
how the PES behaves in the low energy regime, first off 
we should account for the -E-dependent dipole matrix ele- 
ments. Using the Slater-type orbitaO for the correspond- 
ing atomic orbital, we calculate independent dipole ma- 
trix elements. The Slater orbital for R n i(r) is given by 

R nl {r) = (2C)" +1/2 [(2n)!]- 1 / 2 r»- 1 e-^, (16) 

where the orbital exponent is determined by a suit- 
able rule. But in 4/ 1 5rf 1 6s 2 configuration of Ce, the 
Slater orbital for 5d gives actually the poor representa- 
tion compared to a more accurate LSD calculatior£3 for 
the atomic wavefunction for Ce. We adopt therefore the 
same functional form of Eq. (|l6|) , but determine the ex- 
ponent C suitably by comparing with the accurate result, 
i.e. we take Cif = 5-0 and Csd = 2.0 (by a Slater-rule, 
Csd will be 0.75). In principle, the atomic orbital and 
photoelectron basis function having an explicit angular 
momentum channel / should be obtained by solving the 
Schrodinger equation under the same Ce atomic poten- 
tial. But in our discussion the basis function is simply 
assumed to be a spherical Bessel function of I, 

y l E (r) = ^(2E)^MV2Er), (17) 
and its normalization follows 

Jr 2 dr<p l E (r)p E ,(r)=6(E-E>). (18) 
So the dipole matrix elements for 4/-emission is given by 
A f (E) = M f J r 2 drR 4f (r)r^(r), (19) 

and for the 5cZ-conduction band emission, we assume 
A d (e, E) has the simple separable form like 

A d (e,E) = M d A' d (e)A d (E). (20) 
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A^(e) is still assumed to have the shape of V(e) as in the 
last section and A d (E) can be expected to have much 5d 
atomic orbital character if we think of the tight binding 
idea for the corresponding energy band. Thus we assume 
the behavior of A d {E) as 



A d {E) oc j r 2 drR 5d ( 



r)r<p l E 3 (r). 



(21) 



Here it should be noted that the possible Z-channel of 
photoelectrons are I = 2, 4 for 4/-emission and I = 1,3 
for 5d-emission due to the angular momentum selection 
rule, but the major channel will be I — 4 and £ = 3, 
respectively. 




FIG. 3. The behavior of dipole elements - A f (E), A d (E) 
are provided with respect to the photoelectron kinetic energy 
E. Note the different energy scale in two behaviors. The 
absolute values are arbitrary. 
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FIG. 4. The u-dependent valence PES are given for 
4/-emission (upper panel) and 5d-emission (lower panel) . The 
spectra at u) = 0.5 in the upper panel or u) = 4.0 in the lower 
will be so tiny that they are not illustrated. In the respective 
inset, the spectra is normalized to have the same area to give 
the change of shape. The unit of u) is the atomic unit (1 au 
= 27.2 eV). 



In Fig.|3|, we give the dipole matrix behaviors of Af (E) 
and A d (E). This is obtained from very crude calcula- 
tions, but qualitatively quite consistent with Yeh and 
Lindau'sa calculation of photoionization cross sections 
for Ce. We can simply expect from the behaviors of 
dipole elements the general trend is that in the low en- 
ergy, d-emission will be dominant over /-emission, while 
in the high energy, /-emission dominant over (i-emission. 
We first show the calculation result for separate 4/ and 
5d contributions with respect to various photon energies 

Figure [| shows w-dependent 4/- and 5<i-emission, the 
changes of spectral height and shape with u> varied. The 
spectral height will be proportional to the square of 
dipole matrix elements at the corresponding energies and 
the shape related to the behaviors of dipole element. If 
the bd radial wave function does not vary significantly for 
La and Ce, the bottom panel of the figure can be com- 
pared with the PES for Laliil and found to be consistent 
with the experiment. We can also find as w increases 
the spectral shape approaches the sudden approximation 
results (see the insets). 

Now we investigate the total valence PES to which 
both 4/- and 5d-emission contribute with respect to var- 
ious photon energies u>. To parameterize the relative 
strength of two dipole matrix effects, we define Af and 
A d as Af ee A/(£;)|e = 4.o, A d ee max{A d {e, E)}. That 
is, Af is defined as the value of Af(E) at E = 4.0 and 
A d as the value of A d (e, E) at e = 0.0 and E w 0.6 (see 
Fig.|). We give the behaviors of valence PES as w varied 
for \A d /A f \ = 2.0 and |A d /A/| = 6.0, respectively. 
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FIG. 5. The valence PES with respect to photon ener- 
gies are provided. The relative strength of two channels is 
taken as |Ad/A/| = 2.0. Here the solid line is the total 
spectra corresponding to constructive or destructive interfer- 
ence. The dashed line represents 4/-emission, the dotted line 
5d-emission, and the dot-dashed line the interference contri- 
bution. 



In Figs.|| and q, we see that at ui — 0.5 au, the dom- 
inant contributions are from 5d-band emission because 
Af(E) increases slowly compared to A d (e, E), however at 
u> = 4.0 au, most of the spectra in Ce arises from the 4/ 
electrons because A d (e, E) rapidly falls off over E ~ 0.6 
au. Gunnarsson and Sch6nhammerE3 have obtained the 
total emission spectra involving the interference (of 4/ 
and 5d) based on the sudden approximation, but could 
not have discussed these behaviors with respect to w. In 
the experiments, on the other hand, the increasing 4/ 
and decreasing band features witt-, varying u> has been 
used to separate the 4/ structuresEj. That is, in the ex- 
periments, using He resonance lines, to subtract u> = 0.78 
(= 21.2 eV) result from w = 1.5 (= 40.8 eV) result leads 
to approximately a 4/-emission for a moderate value of 
\A d /A f \ (i.e. say for \A d /A f \ ~ 2). 




ENERGY (eV) 



FIG. 6. The valence PES with respect to photon ener- 
gies. The relative strength of two channels is taken as 
| Ad/A/ 1 = 6.0. Notations are same as in FigJ^. 



Below, in Fig.f5| we can see the behaviors of inter- 
ference as D varies. As u) increases, the interference 
become stronger at first and then weaker again. And 
about lo = 1.0 or 2.0 au, we see the strong interfer- 
ence, where the spectra cannot be understood from two 
separate emission spectra. Particularly in the case of 
strong destructive interference, although the 4/-emission 
always comprises two peaks (see the inset of upper panel 
in Fig.Q) , the peak near the Fermi level may be smeared 
until u becomes a bit larger. 




FIG. 7. The contribution of interferences are given with 
respect to the photon energies. In the intermediate energies, 
the interferences are very strong. 



V. LOW ENERGY VALENCE PHOTOEMISSION: 
EFFECTS OF PHOTOELECTRON SCATTERING 

In the low photon energy region, we should consider 
the effects of photoelectron scattering potential induced 
by the hole left by the electron emission as well as the 
dipole matrix behavior. The band emission can also raise 
the shake-up effects like plasmon satellites due to the fluc- 
tuation potentials. Nevertheless, in the present model, 
within 1/Nf expansion, the relevant bases of Eqs.(|^)-([Tc|) 
in the limit of Nf — ► oo do not allow any conduction 
electron-hole excitation. Any shake-up behaviors from 
dielectric responses by band emissions then cannot be 
seen in the taken limit, but in the next higher order of 
1/Nf. For a hole in a localized /-level, however, a small 
number of electrons may undergo measurable shifts in 
response to the potential induced by a holecj. For the 
photoelectron scattering potential, we should go back to 
Eq.([j]) and see the interaction of a /-level impurity elec- 
tron. In this section, we'll confine our discussion only to 
the /-level valence photoemission, i.e. here we do not 
consider the interference with band contributions. In H 
(= Hq + A), it needs noting that the /-electron correla- 
tion is actually a quantity renormalized by the conduc- 
tion electrons, that is, 

u ff n f n f + u fdn f / de<PtuAi> 

= (Uff -U fd )n f nf = Un f n f , (22) 

where we have used rif + ^„ J detl>l u i/j ev is a conserved 
quantity. In the similar way, we can compose the scat- 
tering potential term Vsc which must be added to TLq 

Vsc = V 4f (r)n f + V 5d (r) ]T f de^ ev - V if {r), (23) 

where it should be noted that the initial neutral (ground) 
state is 4/ . Then we have 



G 



V SC = [V Af (r) - V 5d (r)]n f - V if (r) 

= [V 4f (r)-V 5d (r)}(n f ~l)~V 5d (r). 



(24) 



We know V$ d (r) is much broader and weaker than V4 / (r) 
and better to be neglected. So we take Vsc as 



Vsc = V(v)(n f - 1), V(r) = V if {r) - V 5d (v). 



(25) 



Then we cut off the potential by taking the muffin-tin ra- 
dius r m t as the radius of neutral Ce atom, r mt = 3.49 au. 
V4/(r) and V5 d (r) are evaluated from the Slater orbital. 
Now, V(y) is a short range one due to a screening of con- 
duction e?-electrons and have a simple relation between 
V(r) and the intra-atomic Coulomb correlation £/, 



U = drp f (r)V(r) « V(0), 



(26) 



if we assume the /-level charge density p/(r) and /-level 
is quite localized like the core level. Thus V(r) should 
be redefined by ^V(r), where e is a dielectric constant 
chosen to make sure of Eq.(^6|), being due to screening 
by the surrounding, so e is V(Q)/U « 5.24. The behavior 
of V(r) is given in Fig.pl(a). 
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FIG. 8. (a) The photoelectron scattering potential V(r) is 
given as normalized by U. (b) The diagonal and off-diagonal 
parts of scattering potential matrix are given. 

Then we express the scattering potential in terms of 
the photoelectron basis function, 



V sc=Y.J dEdE'V(E,E')^ E ^ E ,, 



(27) 



_ v' 



where the potential matrix element V{E 1 E') are 



V(E,E')= / dv^ Ev {v)V{r)^ E , v (v). 



(28) 



As in calculating the dipole matrix, the photoelectron 
basis function <y9_E^(r) must be obtained by solving the 
Schrodinger equation under the atomic potential. But 
here we use simple spherical Bessel function of / = 4 as 
in the last section. We hopefully expect the essential 
feature will not be spoiled by neglecting the phase shift 
5(k). So the desirable (/^(r) is 

VE u(t) = ^(2S) 1 /4 j4 (y2^ r )y 4m (f.) ; (29) 



and the matrix element V(E, E') is 

V(E, E 1 ) = -{AEE') 1 ^ J drr 2 j 4 (V2Er)V(r)j 4 (V2E J r), 

(30) 



whose explicit behaviors are shown in Fig.[|(b). 

As Vsc added, the total Hamiltonian Tt becomes 
Ho + A + Vsc- Under Tt, the valence PES via /-channel 
can be calculated. The comparison of the results with 
Vsc to those without Vsc (still including the dipole el- 
ements) is provided in Figj9[ The effects of scattering 
potential are quite small as shown in Figj^, which must 
be due to a weak scattering potential in a typical Ce 
compounds. Nevertheless, it's very meaningful to pur- 
sue a general consensus in the valence PES about the 
photoelectron scattering effects. In order to be more in- 
structive, we also investigate the resulting behaviors for 
a slightly different potential whose range is a bit larger 
by 50% (see FigM). 



• = 0.1 




• = 2.0 




1.0 




< = 1.0 




FIG. 9. Effects of scattering potential (solid line) are illus- 
trated by comparing with noninteracting results (dashed line) 
at given photon energies. 
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FIG. 10. Notations are same as in Fig.y. Here the scatter- 
ing potential range has been made a bit larger than in a Ce 
case, r m t = 5.24 au. 



VI. ADIABATIC-SUDDEN TRANSITION FOR 
2-ELECTRON AND N f = 2 

The AIM can be reduced into the two-level model, i.e. 
the whole continuum band is replaced by one level. The 
Hamiltonian H.q we should now consider is 

Ho = e d 1>\j>da + c/ Yl ^f^f" 

<7 <7 

+ v Y,^}^ d - + ^f-) + Un n n fh (3i) 

a 

where nf a = ipt^ipfcr- Then the analogy of the present 
problem with the core level PES for the shake down case 
(having "level crossing" as the hole is created) is more 
evident. The change of relevant electronic levels is given 
schematically before and after the photoemission in the 
following Fig.|i"i"|. 



In both Figs|] and|i"o|, it's notable that there are no ap- 
preciable changes in the peak near the Fermi level, while 
an increase in the ionization peak. This will be under- 
stood from the scattering potential, V(r)(nf — 1). The 
Fermi-level peak is from 4/ 1 — > 4/ 1 , which will not be af- 
fected much by the potential because of rif = 1, but the 
ionization peak is from 4/ 1 — > 4/°. It is important to 
grasp the underlying physics from the spectral changes 
as the potential range is increased from r mt = 3.49 to 
Tmt — 5.24 (1.5 x 3.49). Naturally a longer range poten- 
tial results in more prominent effects in the spectra. In 
two respects, the spectral changes due to the photoelec- 
tron scattering look very similar to the core leveLPES in 
"shake down" case in the previous work of oursEZl. First, 
if we simulate the absorption intensity ratio by the ra- 
tio of the ionization peak to the Fermi-level peak, we 
find the constructive interference between intrinsic and 
extrinsic processes in the low energy regions just as in 
the shake down scenario of core level PES (the ratio is 
increased due to the scattering). The second point is 
the relevant energy scale governing the constructive in- 
terference due to the photoelectron scattering. In FigJ|, 
we see the maximal scattering effects around lu ~ 2.0, 
while, in Fig.[l0], the energy scale giving the maximal ef- 
fects is around uj ~ 1.0. That is, the governing energy 
scale is decreased by an increased range. In the core level 
PES, we found the relevant energy (E) can be directly 
related by the potential range (R) as E = l/(2i? 2 ), which 
tempts the application of the criteria to the present sys- 
tem. And it then can be a natural motivation that we 
make an analysis for the adiabatic-sudden transition in 
this system and try to answer if the criteria found previ- 
ously in core level PES can be still valid in this valence 
PES or not. This question is extensively discussed in the 
next section, where we propose the simplified two-level 
model for the sake of simplicity. 



r 



p 



before 



after 



FIG. 11. Schematic view of the relevant configurations in 
initial and final stage. Here we assume 2e/ + U w 0. Note 
there's a level crossing before and after the emission. 

We can introduce three states |/°), If 1 ), and \f 2 ) as 
follows, 

l/°)=^W||0). 



I/ 1 ) = ^[^/;I°)-^} T I 
|/ 2 )=VV}J0), 

then we can express TIq in these bases, 

V 
H.: | I Ai V i + •!<„ 
V 2Ae + U 



(32) 



where Ae = ej — €d and V — \/2V . For simplicity, we 
will put 2Ae + U = 0. The ground state corresponding 
to the initial state of photoemission, I^q) lS 



l*o> 



V 



im^imi/ 2 ) 



(33) 



where A = iAe — iy Ae 2 + 8V 2 and its energy E is 



E = A + 2e d . 



(34) 
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The final states of the target are given by the following 
another set of bases, 

|/°;<7}=4 CT |0), 
|/V>=^|0>, 

then the Hamiltonian H.q with one /-electron emitted is 



Ho 



( V/y/2\ 
\V/y/2 Ae J 



(35) 



Note ^j fa \f°) = 0, V/.l/ 1 ) = ^<7\f°;-a), and 

V'/o-l/ 2 ) = o - !/ 1 ;— cr) and a is ±1. Here a is actually 
a redundant parameter. The possible final target states 
will be given by the eigenstates of 7Yq m Eq. ( p5| ) , 



|*i;<t) = cosv3|/°;o-) - sin^l/ 1 ; ct), 



1*2; <t) = sin^|/°;a) +cos(^|/ 1 ;a), 



(36) 



(37) 



where Ei = ±Ae =F ^VAe 2 + 2V 2 + e d (SE = 

V Ae 2 + 2V 2 ), and the parameter ip (-| < ip < -|) is 
determined by 



1 . / — Ae 

tamp = — (vr + 2 — u>), w = 
V 2 



(38) 



Also noticeable is 



il> f< r\*o) = -7= [-sin0|/°;-a) +cos0|/ 1 ;- ( r>J , (39) 



where the parametric angle 9 (f < < |) is 

cot0 = -^p(Vw 2 + 8 + w). 
2y 2 



(40) 



Now we consider the optical activation Hamiltonian A, 
A = Y, Mk1>t,1>f<r, (or A = J dEA{E)^ fa \ . 

(41) 



Within the first order perturbation theory, the photoc- 
mission matrix element M(i, ka) (i = 1, 2) will be 



M(i, ka) 
(^i;-a\ip k c 



1 + V 



1 



iM k + 2J c ij m j 



E-Ho-T + iri 
V w M k 



A|* ) 



■? u + E -Ej 



1 I./2 



j fc' ~ 1 ^ 2 

where the scattering potential Vsc is taken as 



IT] 



(42) 



y 5 c - V(r)(n/ - 1), V kk , = J dr<p%„{*)V{r)vv a {r) 

(43) 

and thus 

mi = (#i;-oi^ /ff |#o), = -a\n f \^f, -a) - 
That is, the coefficients are 

mi = — -= sin(</? + 8)a, m 2 = — ;= cos(ip + 8)a, 
V 2 v 2 



Cll 



COS C22 



sin (p, C12 = c 2 i 



sin tp cos (/p. 



If we consider a ratio between the main and the satel- 
lite absorption intensity divided by a nonintcracting case 
r(aj)/r (uj), r(w)/r Q (w) is 



l + sin 2 ^ (I) 



(44) 



:.n 2tp sin((p+£ 
2cos(ip+0) 



) V 171 / li+iEA 



COS 2 V^Ffc 



/ d>+<5_E \ 
1 \ E ) 



sin 2tp cos(tp+#) y iTi 

2sin( V +0) £ ^ \ B 



where fci = J2(w + Eq - E 1) , Cu = E 2 — E a (threshold 

energy for the satellite), and we have extended the model 
matrix elements M k and Vkk' used in the previous core 
level case, 



E 



Vkk>M k . 
e - ey + it] 



1 f c 
Fk(e) = - / 



ZM k F k {e/E), 
E 



x 10 dx 



(45) 



[1 + x 5 ] 2 l + ^fc-a;) 2 [x 2 -e-ir]} 

(46) 



M fc = 



(Rk) 5 



1 + (i?fc) 5 ' 



(47) 



{R 2 kk'f 



R [1 + (i?fc) 5 ] [1 + (i?fc') 5 ] [1 + R 2 (k- fc') 2 



(48) 



where i? is the characteristic length scale of the system 
directly related to the potential range and E = 1/2R 2 . 
Here it is found that from Eq.(Q), r(w)/r (w) can be 
written essentially in the same mathematics as in the 
core level case. In the following Fig.[l^, we give the be- 
haviors of F(e). 
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FIG. 12. The function F(e k ) denned in Eq.©. Both the 
real and imaginary parts are given. 

Similarly to the core case, we always have an over- 
shoot behavior in r(u>)/ro(uj) in the low energy limit, 
when SE = 0, 



ro(u>th) 



rp(rt\ V sin yj sin 9 
r y V )~E cos(y+e) 
pfr\\ V cos ip sin 6 



> 1, 



(49) 



where F(0) = 0.052286 and f < ^ < f , ?< < f , 
and ? < 9? + < n should be noted. In Fig.|13|, we show 
r(u)/ro(itj) as a function of ui/E for a few values of V/l? 
(u> = u>- UJ t h)- 




0.2 



E = 0.2 (V/B = 1.0) 
£ = 0.3 (V/E = 0.7 
£ = 0.4 (V"/£ = 0.5) 
E = 0.5 (V/B = 0.4) 
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j/E 



15 



20 



FIG. 13. The ratio r(oj)/ro((jj) as a function of uj/E for 
several values of V/E. Ae/V = —2.0 is taken. 

Fig.|l3| shows as ui increases the ratio also increases and 
reachee a maximum. From the arguments in our previous 
workO, we see roughly R~R /3 and V ~ 3V(0)/2 and 
for Ce, R ~ 1 and U - 0.3, which leads to VyjS - 0.6. 
Most notable is that the curves have the universal fea- 
ture independent of V/E, i.e. the overshoots disappear 
at about uj/E ~ 10 in all relevant parameter region. This 
means the adiabatic-sudden transition depends only on 
E, that is, R, even if the amplitude of overshoot is from 
V/E. Beyond the first order perturbation, the overshoot 
range will be reduced due to the multiple scattering, but 
the universal behavior not changed. This coiiclusion is 
exactly identical to that in the core level casal3 and im- 



plies the same criteria can be applied also to the valence 
PES case. 



VII. CONCLUSION 



We have studied the valence photoemission spectra in 
the Anderson impurity model aiming at Ce compounds. 
Using the time-dependent formulation and l/Nf expan- 
sion, we can treat the problem exactly up to 0(1/Nf) . 
For Ce compounds, Nf — oo can be a good approxi- 
mation for Nf = 14. Within the formalism, to evaluate 
the photoemission spectra is corresponding to solving the 
time-dependent Schrodinger equation. 

To investigate the low energy photoemission spectra, 
we should consider the dipole matrix and photoelectron 
scattering matrix additionally compared to the sudden 
approximation valid in high energy limit. In view of ex- 
periment, the valence PES always consist of /-emission 
and band emission. So we considered both dipole matrix 
elements having explicit E'-dependencies and obtained 
the total spectra as well as two separate spectra with 
respect to the photon energies. The relative strength 
and sign of two dipole elements are crucial in the total 
spectra. Due to differences in the energy scales of 4/ 
and 5c? dipole elements, the general trends are that, in a 
very low energy (w < 1.0 au), the 5d emission is domi- 
nant, while the 4/ emission increases and dominates over 
the 5d emission in a high energy (u> > 2.0 au). It is also 
found that the interference effects of /-emission and band 
emission (from d) are highly nontrivial especially in the 
intermediate energy region. In the intermediate region 
(uj ~ 1.0 or 2.0 au depending on the relative strength), 
due to a strong peak of the interferences near the Fermi 
level, the Kondo resonance peak of 4/ emission may be 
smeared in the case of destructive interference. The con- 
structive or destructive interference will be determined 
by the relative sign. 

We also studied the effects of scattering potential for 
the /-electron emission. Scattering effects on the band 
emission is neglected in the infinite Nf limit. The poten- 
tial matrix also includes the kinetic energy dependencies. 
It's a general result that, in the case where the photo- 
electron couples to the localized excitation, the arrival 
at the sudden transition is much faster compared to the 
extended excitation. Effects of scattering potential gives 
the similar spectral changes to the core level PES in shake 
down case if we assign the Kondo resonance and ioniza- 
tion peak to the main and satellite peak, respectively. 
Behaviors of two-peak ratio is reminiscent of the previous 
analysis of core level PES. Therefore, this can be a moti- 
vation to do a further analysis for the adiabatic-sudden 
behavior. We can then cast a question, "The criteria 
found in core level case can be also valid in the valence 
case?" To explore it, first we simplify the model into just 
a two-level one, where the whole conduction band is re- 
placed by one level. Then it is found the intensity ratio 
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of two peaks can be written actually in the same way 
as in the core case. Through the same analysis, we can 
find that also in the valence PES, the sudden transition 
happens on the energy scale of E = 1 / (2i? 2 ) , where R is 
a typical length scale of the scattering potential. 

Finally, we would like to make it clear the limit or 
shortcomings of our present model. Within the model, 
we cannot yet arrive at the realistic photoemission in 
a true solid. The AIM includes the continuum band de- 
scribing a crystal, which has a surface. Then the position 
of a surface with respect to the created hole will enter as 
an important parameter. In a true solid, we never reach a 
vanishing extrinsic scattering because the corresponding 
cross section goes to zero, while the photoelectron is emit- 
ted from an infinite depth. And to correctly predict the 
amount of photocurrent, we also should use the damped 
wavefunctions, not undamped partial wave. These addi- 
tional extrinsic effects may be as important as the hole- 
induced scattering potential treated here. Therefore, it 
should be understood that the scope of this paper is to 
explore the effects of hole- induced scattering potential in 
the valence PES in AIM describing the Ce compounds 
apart from additional ingredients. 
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